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Abstract. I shortly review the present status of the theoretical calculations of e'/e and the compar- 
f~>) ■ ison with the present experimental results. I discuss the role of higher order chiral corrections and 

| in general of non-factorizable contributions for the explanation of the AI = 1/2 selection rule and 

direct CP violation in kaon decays. Still lacking reliable lattice calculations, analytic methods and 
phenomenological approaches are helpful in understanding correlations among theoretical effects 
and the experimental data. Substantial progress from lattice QCD is expected in the coming years. 
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INTRODUCTION 

> 

The results obtained in the last two years by the NA48 [1] and the KTeV [2] collabora- 
tions have marked a great experimental achievement, establishing some 35 years after 
the discovery of CP violation in the neutral kaon system [3] the existence of a much 
smaller violation acting directly in the decays: 

O 

a: 



Ref c'/^- J (15.3±2.6)xl0- 4 (NA48) 

lt/tJ ~ \ (20.7 ±2.8) x 10- 4 (KTeV). UJ 

The average of these results with the previous measurements by the NA31 collaboration 
at CERN and by the E731 experiment at Fermilab gives 

^- Re(e'/e) = (17.2 ±1.8) x 1(T 4 (2) 

While the Standard Model (SM) of strong and electroweak interactions provides an 
economical and elegant understanding of indirect (£) and direct (£ ; ) CP violation in term 
of a single phase, the detailed calculation of the size of these effects implies mastering 
strong interactions at a scale where perturbative methods break down. In addition, direct 
CP violation in K — > nn decays arises from a detailed balance of two competing sets 
of contributions, which may hopelessly inflate the uncertainties related to the relevant 
hadronic matrix elements in the final outcome. All that makes predicting e'/e a complex 
and challenging task [4] . 

Just from the onset of the calculation the presence in the definition of e'/e , written as 



e'_ 1 f ({%%) I=2 \Xw\K L ) ((rni)i= 2 \Xw\Ks) \ 
e V2\((^)i=o\^w\K L } ((Ttn)i=o\rtw\Ks) J 



(3) 



of given ratios of isospin amplitudes warns us of a longstanding and still unsolved 
theoretical "problem": the explanation of the AI = 1/2 selection rule. 



The AI = 1/2 selection rule in K — > tttt decays is known since 45 years [5] and it 
states the experimental evidence that kaons are 400 times more likely to decay in the 
1 = two-pion state than in the / = 2 component (u) = A2/A0 ~ 1/22). This rule is not 
justified by any symmetry argument and, although it is common understanding that its 
explanation must be rooted in the dynamics of strong interactions, there is no up to date 
derivation of this effect from first principle QCD. 

Given the possibility that common systematic uncertainties may a-priori affect the 
calculation of e'/e and the AI = 1/2 rule (see for instance the present difficulties in 
calculating on the lattice the "penguin contractions" for CP violating as well as for CP 
conserving amplitudes [6]) a convincing calculation of e'/e must involve at the same 
time a reliable explanation of the AI = 1/2 selection rule. Both observables indicate the 
need of large corrections to factorization in the evaluation of the four-quark hadronic 
transitions. Among these corrections Final State Interactions (FSI) play a substantial 
role. However, FSI alone are not enough to account for the large ratio of the / = 
over 1 = 2 amplitudes. Other sources of large non-factorizable corrections are therefore 
needed for the CP conserving amplitudes, which might affect the determination of 
e'/e as well. As a consequence, a self-contained calculation of e'/e should also address 
the determination of the CP conserving K ^ nn amplitudes. 



The Operator Product Expansion (OPE) provides us with a very effective way to address 
the calculation of hadronic transitions in gauge theories. The integration of the "heavy" 
gauge and matter fields allows us to write the relevant amplitudes in terms of the 
hadronic matrix elements of effective quark operators and of the corresponding Wilson 
coefficients (at a scale fi), which encode the information about those dynamical degrees 
of freedom which are heavier than the chosen renormalization scale. With the SM flavor 
structure AS = 1 transitions are effectively described by 



The entries V\j of the 3 x 3 Cabibbo-Kobayashi-Maskawa (CKM) matrix describe the 
flavour mixing in the SM and x = —V t dV t */V u dV* s . For fi < m c (q = w, d, s), the relevant 
quark operators are: 



OPE: AN "EFFECTIVE" APPROACH 



XAs =i = -jfudV^Y^) +iyi<ji)\Qi{n) ■ 



(4) 
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Current-Current 



23,5 
24,6 




Gluon "penguins" 



(5) 



27,9 
28,10 




Electroweak "penguins 



Current-current operators are induced by tree-level W-exchange whereas the so-called 
penguin (and "box") diagrams are generated via an electroweak loop. Only the latter 
"feel" all three quark families via the virtual quark exchange and are therefore sensitive 
to the weak CP phase. Current-current operators control instead the CP conserving 
transitions. This fact suggests already that the connection between e'/e and the AI = 1/2 
rule is by no means a straightforward one. 

Using the effective AS = 1 quark Hamiltonian we can write e'/e as 



e' m G F co 



g' ^lTrT, Im ^ 

2|£|ReA 



n --n 2 

(0 



(6) 



where 

n = ^L^Re(a) (i-a IB ) 
n 2 = iL^Re<a> 2 , 

and (<2;) = (%%\Qi\K). The rescattering phases 80,2 can be extracted from elastic it — it 
scattering data [7] and are such that cos 80 ^ 0.8 and cos 82 ^ 1. Given that the phase 
of £ (0 e ) is approximately n/4, as well as 80 — 82, § = § + 82 — 80 — e turns out 
to be consistent with zero. While Gp, (Q, |e| and ReAo are precisely determined by 
experimental data, the first source of uncertainty that we encounter in eq. 6 is the value 
of lmX t = lm(V t * s V t d), the combination of CKM elements which measures CP violation 
in AS = 1 transitions. The determination of lmX t depends on B-physics constraints 
and on £ [8]. In turn, the fit of £ depends on the theoretical determination of Bk, the 
K° — K° hadronic parameter, which should be self-consistently determined within every 
analysis. The theoretical uncertainty on Bk was in the past the main component of the 
final uncertainty on lmX t . The improved determination of the unitarity triangle coming 
from B-factories and hadronic colliders [9] has weakened and will eventually lift the 
dependence of ImX; on B K , allowing for an experimental measurement of the latter 
from £ . Within kaon physics, the decay Kl — > 7t°VV gives the cleanest "theoretical" 
determination of lmX t , albeit representing a great experimental challenge. At present, a 
typical range of values for lmX t is (0.94 - 1.60) x 10~ 4 [10]. 

We come now to the quantities in the square brackets. While the calculation of 
the Wilson coefficients is well under control, thanks primarily to the work done in 
the early nineties by the Munich [11] and Rome [12] groups, the evaluation of the 
"long-distance" factors in eq. 7 is the crucial issue for the ongoing calculations. The 
isospin breaking (IB) parameter Qib, gives at the leading-order (LO) in the chiral 
expansion a positive correction to the A2 amplitude (proportional to Aq via the n° — T] 
mixing) of about 0.13 [13]. At the next-to-leading order (NLO) the full inclusion of the 
7t° - T| - t\' mixing lift the value of to 0.16 ± 0.03 [14]. On the other hand, the 
complete NLO calculation of IB effects beyond the n° — r\ — r\' mixing (of strong and 
electromagnetic origin, among which the presence of AI = 5/2 transitions) involves a 
number of unknown NLO chiral couplings and is presently quite uncertain. Dimensional 
estimates show that IB effects may be large and affect e'/e sizeably in both directions 
[15]. Although a partial cancellation of the indirect (Aro) and direct AQib NLO isospin 
breaking corrections in eq. 6 may reduce their final numerical impact on e'/e , we must 
await for further analyses in order to confidently assess their relevance. At present one 
may use Qib = 0.10 ±0.20 [16, 17, 18] as a conservative estimate of the IB effects. 



FIGURE 1. Anatomy of e'/e in the Vacuum Saturation Approximation. In light (dark) gray the positive 
(negative) contributions of the effective four-quark operators are shown with proportional weight. 

The final basic ingredient for the calculation of e'/e is the evaluation of the K — > 7UT. 
hadronic matrix elements of the quark operators in eq. 5. A simple albeit naive approach 
to the problem is the Vacuum Saturation Approximation (VSA), which is based on two 
drastic assumptions: the factorization of the four quark operators in products of currents 
and densities and the saturation of the intermediate states by the vacuum state. As an 
example: 

{% + lC\Q 6 \K Q ) = 2{%-\uy 5 d\Q){%^\MK )~2{%^%-\dd\Q){0\^ 5 d\K Q ) 

+ 2 [(0\ss\0) - (0\dd\0)] (7r + ^|5y 5 J|^°) (8) 

The VSA does not exhibit a consistent matching of the renormalization scale and 
scheme dependences of the Wilson coefficients and it carries potentially large systematic 
uncertainties [4] . On the other hand it provides useful insights on the main features of the 
problem. A pictorial summary of the relative weights of the contributions of the various 
operators to e'/e , as obtained in the VSA, is shown in Fig. 1. 

As we have already mentioned, CP violation involves loop-induced operators (Q$ — 
Qio). From Fig. 1 one clearly notices the potentially large cancellation among the strong 
and electroweak sectors and the leading role played by the gluonic penguin operator Q(, 
and the electroweak operator Q%. Tipical range of values for e'/e , obtained using the 
VSA, are shown in Fig. 2 together with the three most updated predictions available 
before 1999 (when the first KTeV and NA48 results became known) [19, 20, 26]. The 
fact that the cancellation among the strong and electroweak sectors turns out to be quite 
effective (at least in the VSA) warns us about the possibility that the uncertainties in 
the determination of the relevant hadronic matrix elements may be largely amplified in 
the calculation of e'/e . It is therefore important to asses carefully the approximations 
related to the various parts of the calculations. In particular, the analysis of the problem 
suggests that factorization may be highly unreliable. 



4 





3 


CO 


2 


O 








X 






1 













-1 



NA31 



E731 



1 



VSA 
HVNDR 



Muenchen 
(1996) 



Roma 
(1997) 



Trieste 
(1997) 



PRE-DICTIONS 

FIGURE 2. The l-o results of the NA31 and E731 Collaborations (early 90's) are shown by the 
gray horizontal bands. The old Miinchen, Roma and Trieste theoretical predictions for e'/e are depicted 
by the vertical bars with their central values. For comparison, the VSA estimate is shown using two 
renormalization schemes. 



BEYOND FACTORIZATION 

The dark gray bars in Fig. 2 represent the results of three calculations of e'/e which are 
representative of approaches that in principle allow us to go beyond naive factorization. 
They are based from left to right on the large N c expansion [19, 21], on lattice regu- 
larization [20, 22], and on phenomenological modelling of low-energy QCD (the chiral 
quark model) [23, 24, 25, 26]. 

The experimental and theoretical scenarios have changed substantially after the first 
KTeV data and the subsequent NA48 results. Fig. 3 shows the present experimental 
world average for e'/e compared with the revised or new theoretical calculations that 
appeared during the last year. Without entering into the details of the results (for a 
short summary see [27]) they all represent attempts to incorporate non-perturbative 
information into the calculation of the hadronic matrix elemnts, whether it is based 
on the large N c expansion (Miinchen [28], Dortmund [29], Beijing [30], Taipei [31], 
Valencia [32]), phenomenological modelling of low-energy QCD (Dubna [33], Tri- 
este [34, 35], Lund [36]), QCD Sum Rules (Montpellier [37]) or, finally, lattice regu- 
larization (Roma [6]). 

Overall most of the theoretical calculations are consistent with a non-vanishing pos- 
itive effect in the SM (with the exception of some very recent lattice results on which I 
will comment shortly). 

At a closer look however, if we focus our attention on the central values, many of the 
predictions prefer the 10~ 4 regime, whereas only a few of them stand above 10~ 3 . Is this 
just "noise" in the theoretical calculations? Actually, without entering the many details 
on which the estimates are based, most of the aforementioned difference is explained in 
terms of a single effect: the different size of the hadronic matrix element of the gluonic 
penguin Q(, obtained in the various approaches. In turn this can be understood in terms 
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FIGURE 3. Recent theoretical calculations of e'/e are compared with the combined l-o average of the 
NA31, E731, KTeV and NA48 results (e'/e = 17.2 ± 1.8 x 1(T 4 ), depicted by the gray horizontal band. 

of sizeable higher order chiral contributions (NLO in the l/N c expansion) to the 1 = 
amplitudes. 

This effect was stigmatized before the latest experimental round by the work of the 
Trieste group [25, 26], and appears clearly in the comparison of the leading l/N c and 
lattice results with the chiral quark model analysis in Fig. 2. The chiral quark model 
approach, together with the fit of the CP conserving amplitudes to phenomenologically 
normalize matching and model parameters, allows one to carry a complete calculation of 
the hadronic matrix elements beyond the leading order in the chiral expansion (including 
the needed local counterterms). Non-factorizable chiral contributions (not included in 
leading l/N c or lattice calculations) were shown to produce a substantial enhancement 
of the 1 = transitions lifting the expectation of e'/e at the 10~ 3 level. 

Since then a number of groups have attempted to improve the calculation of K — > 7UT. 
matrix elements in a model independent way. The present status of the art is summarized 
in Table 1. Due to the leading role played by and Q% we may write a simplified 
version of eq. 6, namely 

which in spite of the fact that it should "not be used for any serious analysis" [28] gives 
an effective and practical way to test and compare different calculations. 

The B-factors B\ = (Qi) / (<2/)vsa represent a convenient parametrization of the 
hadronic matrix elements, albeit tricky, in that their values are in general scale and 
renormalization- scheme dependent, and a spurious dependence on the quark masses is 
introduced in the result whenever quark densities are involved. The latter is the case 
for the Q(, and gs penguins. As a consequence the VSA normalization may vary from 
author to author thus introducing systematic ambiguities. By defining the VSA matrix 



TABLE 1. Comparison of various calculations of penguin matr ix ele- 
ments. The data marked by the star are rescaled by a factor ^/3/2, to ac- 
count for a different definition of the isospin matrix elements. 



Method B 6 (NDR) Z?g (NDR) 



Lattice [CP-PACS] 


<0.3 


-0.9 


[38] 


Lattice [RBC] 


-0.4 


- 1 


[39] 


Lattice [APE] 




0.58 ±0.06* 


[40] 


Lattice [SPQcdR] 




0.56 ±0.07* 


[41] 


Large N C +LMD (x-limit) 




4.0±1.2* 


[42] 


Dispersive+data (x-limit) 




2.5 ±0.8 * 


[43] 


Dispersive+data (x-limit) 




1.4 ±0.6 


[44] 


Large N c + data 


1.0 ±0.3 


0.8 ±0.2 


[28] 


NLO 1 /N c CHPT 


1.5 ~ 1.7 


0.4-0.7 


[29] 


NLO 1 /N c ENJL (X-limit) 


2.9 ±0.5 


1.5 ±0.2 


[36] 


NLO x-Quark Model 


1.5 ±0.4 


0.84 ±0.04 


[26, 34] 


Large AT e + FSI 


1.55±0.10 


0.92 ±0.03 


[32] 



elements at the scale fi = 2 GeV we obtain [4] 

((™)i=2\Q&\K°) YSA = V6fm 4 K (m s + m d )- 2 ~ 1.1 GeV 3 , 

<7t7t| Q 6 1^°> VSA / < (7t7T:) /= 2| 08 1^°> VS A = (/*"/*)//* - "0.63 , 

where I have used (m s + ntd) (2 GeV) = 110 MeV and the chiral value / = 86 MeV for 
the octet decay constant. 

(2) 

It is known that 5g and B\ are perturbatively very weakly dependent on the renor- 
malization scale [4]. Therefore it makes sense to compare the B's obtained in different 
approaches, where the matrix elements (Qi) are computed at different scales. In Table 1 
I collected some recent results for the relevant penguin matrix elements coming from 
various approaches, having care of normalizing the results in a homogeneous way (as 
far as detailed information on definitions and renormalization schemes was available). 

As a guiding information, taking ImA* = 1.3 x 10~ 4 , the present experimental central 

value of e'/e is reproduced by (1 — Q.m) — 0.4 B^ pa 1. 

The most important fact is the first evidence of a signal in lattice calculations of 
(%%\Q 6 \K), obtained by the CP-PACS [38] and RBC [39] collaborations. Both groups 
use the Domain Wall Fermion approach which allows to control the chiral symmetry 
on the lattice as a volume effect in a fifth dimension. This approach softens in principle 
the problem of large power subtractions which affects the lattice extraction of / = 
amplitudes (penguin contractions). Still the (n\Qj\K) transition is directly computed on 
the lattice and chiral perturbation is needed to extrapolate to the physical amplitude. The 
two groups obtain comparable values of the and Q% matrix elements leading both to 
a negative e'/e . On the other hand the calculations are at an early stage and do not fully 
include higher order chiral contributions which may be responsible for the dynamical 
enhancement of the gluonic penguin, as large N c approaches beyond LO and the Chiral 
Quark Model suggest. 

Important results have been obtained using data on spectral functions in connection 
with QCD sum rules and dispersive relations in the attempt to obtain model independent 



information on the relevant matrix elements. These approaches have produced as of 
today calculations of Q% (in the chiral limit) which are substantially larger than the 
factorization (and lattice) results. In addition there is still substantial disagreement 
among different methods, in particular between the result of ref. [42], which is based 
on large N c and Lowest Meson Dominance (LMD), and that of ref. [44] obtained via 
dispersive relations and data. 

Calculations which sofar have gone beyond LO in the large N c expansion, based on 
chiral perturbation theory and/or models of low-energy QCD, have shown the crucial 
role of higher order non-factorizable corrections in the enhancement of the 1 = matrix 
elements [26, 29, 32, 36]. Chiral loop corrections drive the final value of e'/e in the 
ballpark of the present data. However the calculation of higher order chiral corrections 
cannot be fully accomplished in a model independent way due to the many unknown 
NLO local couplings. In the chiral quark model approach all needed local coupling are 
computed in terms of quark masses, meson decay constants and a few non-perturbative 
parameters as quark and gluon condensates. The latter are determined self-consistently 
in a phenomenological way via the fit of the CP conserving K nn amplitudes, thus 
encoding the AI = 1/2 rule in the calculation [25, 34]. The analysis shows that the role 
of local counterterms is subleading to the chiral logs when using the Modified Minimal 
Subtraction as opposed to the commonly used Gasser-Leutwyler prescription. 

Among non-factorizable chiral corrections FSI play a leading role. As a matter of 
fact, one should in general expect an enhancement of e'/e with respect to the naive 
VSA due to FSI. As Fermi first argued [45], in potential scattering the isospin 1 = 
two-body states feel an attractive interaction, of a sign opposite to that of the / = 2 
components thus affecting the size of the corresponding amplitudes. This feature is 
at the root of the enhancement of the 1 = amplitude over the 1 = 2 one and of the 
corresponding enhancement of e'/e beyond factorization. An attempt to resum these 
effects in a model independent way has been worked out by the authors of ref. [32], 
using a dispersive approach a la Omnes-Mushkelishvili [46]. Their analysis shows that 
resummation does not substantially modify the one-loop perturbative result and, as 
it appears from Table 1, a 50% enhancement of the gluonic penguin matrix element 
is found over the factorized result. However, the analysis suffers from a sistematic 
uncertainty due to the value of the factorized amplitude at the chosen subtraction point 
which is identified with the large N c result [47] . Even when the authors in the most recent 
work match the dispersive resummation with the perturbative one-loop calculation, thus 
including l/N c effects, again a systematic uncertainty remains hidden in the unknown 
polinomial parts of the local chiral counterterms. In this respect it is important to stress 
that gauging the scheme dependence of the result by varying the renormalization scale in 
the chiral logs cannot account for the arbitrariness related to the chiral-loop subtraction 
(a model-independent complete calculation is missing). 

Finally, it has been recently emphasized [48] that cut-off based approaches should 
pay attention to higher-dimension operators which become relevant for matching scales 
below 2 GeV and may represent one of the largest sources of uncertainty in present 
calculations. The results of refs. [43, 44] include these effects. The calculations based 
on dimensional regularization may be safe if phenomenological input is used in order to 
encode in the relevant hadronic matrix elements the physics at all scales (this is done in 
the Trieste approach). 



In summary, while some model dependent calculations suggest no conflict between 
theory and experiment for e'/e , a precise and "pristine" prediction of the observable is 
still quite ahead of us. 



OUTLOOK AND CONCLUSIONS 

Higher-order chiral corrections are taking the stage of K — > %% physics. They are needed 
in order to asses the size of crucial parameters (as Qib) and the effect of non-factorizable 
contributions in the penguin matrix elements. 

Lattice, as a regularization of QCD, is the first-principle approach to the problem. 
However, lattice calculations still heavily depend on chiral perturbation theory [49]. 
Presently, very promising developments are being undertaken to circumvemt the techni- 
cal and conceptual shortcomings related to the calculation of weak matrix elements (for 
a recent survey see ref. [50]). Among those are the Domain Wall Fermion approach [51] 
which allows us to decouple the chiral symmetry from the continuum limit, and the very 
interesting observation that the Maiani-Testa theorem [52] can be overcomed using the 
fact that lattice calculations are performed in finite volume [53], thus allowing for the 
direct calculation of the physical K — > nn amplitude on the lattice. All these develop- 
ments need a tremendous effort in machine power and in devising faster algorithms. 
Preliminary results for lattice calculations of both e'/e and the Al = 1/2 selection rule 
are already available and others are currently under way [41]. 

In the meantime analytical and semi-phenomenological approaches have been cru- 
cially helpful in driving the attention of the community on some systematic short- 
comings of "first-principle" calculations. The amount of theoretical work triggered by 
the NA48 and KTeV data promises rewarding and perhaps exciting results within the 
next few years. 
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